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Abstract 

We show that the formally exact expression for the free energy (with a non- 
relativistic Hamiltonian) for the correlated metal generates the Poisson equa- 
tion within the saddle-point approximation for the electric potential, where 
the charge density automatically includes correlations. In this approximation 
the problem is reduced to the self-consistent periodical Anderson model (SC- 
PAM). The parameter of the mixing interaction in this formulation have to 
be found self-consistently together with the correlated charge density. The 
factors, calculated by Irkhin, for the mixing interaction, which reflect the 
structure of the many-electron states of the /-ion involved, arise automati- 
cally in this formulation and are quite sensitive to the specific element we are 
interested in. We also discuss the definitions of the mixing interaction for the 

mapping from ab initio to model calculations. 
71.10.-w,71.27.+a,71.28.+d,71.20.Eh 
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I. INTRODUCTION 



The local density approximation (LDA) to density functional theory (DFT) is working 
surprisingly well in many cases for which is is expected not to work at all. One of its 
most important features is that the self-consistent procedure provides a quite accurate dis- 
tribution of the charge density, calculated from the Kohn-Sham equationB which is more 
accurate than the Poisson equation It is also important to understand why the form of the 
potential which has been derived from the theory of a homogeneous electron gas with the 
charge density which includes contributions from the localized electrons (the extreme case 
of non-homogeneity) works so well. Below we will show that, at least, the Poisson equa- 
tion can be formulated also in the case of strongly correlated system, where some part of 
the electrons are either fully or partly localized. One the one hand, physically it is clear 
that the role played by the localized (core) electrons in the formation of the potential is 
the screening of the nuclear potential experienced by the conduction electrons (which con- 
tribute to the cohesive energy). The localized electrons do not contribute to the cohesive 
energy and, therefore, the error, coming from an insufficient description of the core electrons 
is not so essential for the properties derived from a calculation of the total energy at zero 
temperature. On the other hand, the experience accumulated using the LDA-DFT shows 
that it fails to describe the properties, which require information about (quasi) localized 
electrons, like photo-electron spectroscopy experiments, or exchange interaction between lo- 
calized electrons in magnetic insulators, semiconductors, etc. In these cases, methods either 
completely based on the field theory or the ones combining the field-theory and DFT, are 
desirable. The models often used are the Anderson impurity and periodical models and here 
we will discuss mainly the periodical model. There are two ways for providing methods 
combining DFT and field theory. The first way is the following. The field operators can be 
constructed using the functions generated within a LDA-DFT calculation and a correction, 
constructed from the difference TC mt = T[f uU — 7i LDA , can be used to correct the electronic 
structure generated by the initial LDA calculation!!. The motivation for using the LDA-DFT 
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basis here is that one can expect compensation by the LDA potential from those part of 
the self-energy for the conduction electrons, which corresponds to a static random phase 
approximation!. This way, however, involves complex calculations of the matrix elements 
of the Coulomb interaction and frequency-dependent magnitudes, like the self-energy, that 
in practice such calculations are quite hard to performB!. The other way, much more often 
used, is some mapping to the Anderson or Hubbard models. Then, two difficulties arise. 
First, the question about double counting of some of the interactions, and also how to calcu- 
late the parameters of the model which is chosen for treating correlations. Second, the model 
calculations often involve a redistribution of spectral weights between low- and high-energy 
regions and a redistribution of the charge density caused by it. The latter is never taken 
into account in the model calculations. It is especially important since any redistribution 
of charge involves a large Coulomb energy. This is especially important when the Anderson 
model is used for discussing magnetic properties. These properties are determined by the 
effective exchange integral ~ V 2 /U (where V is the mixing and U is the Hubbard on-site re- 
pulsion) and, therefore, involves small energies. Calculations of this small energy difference, 
having neglected a possibly greater Coulomb energy, can easily lead to a misinterpretation 
of the experimental data. It is also important that the mixing interaction is representation 
dependent and, therefore, for the description of a real system within a model it is important 
to define clearly what is mixing interaction for that special case. This shows the need for 
a formulation which allows for a self-consistent calculation of the parameters of the model 
together with the charge density. Such an opportunity arise in a natural way if one starts 
with the full Hamiltonian and treat the single-ion Coulomb interactions in some approxima- 
tion which takes into account the strong local electron correlations. Such a scheme has been 
suggested earlier!, however, the strong electron correlations (SEC) have been treated within 
the slave-boson technique which at present seems to be unsatisfactory for reasons which we 
will not discuss here (see, Reffl). Here we will discuss the parameters of the Audersou model 
within the same scheme, using the definitions of the operators in a non-orthogonal basis set 
different from Ref.!, but coinciding with the ones used in Ref.i and!. We have discussed a 
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way to calculate the single-site Coulomb repulsion parameter, Hubbard U, earlier^. Here 
we pay attention mainly to the mixing interaction. The organization of the paper is as 
follows. In section II we rewrite the Hamiltonian in an non-orthogonal basis set, and con- 
struct the many-electron operators. Using a saddle point approximation we get the Poisson 
equation for the SEC system. In section III we show how the periodical Hubbard- Anderson 
model appears using the results of section II. In section IV we discuss the mixing parameter 
entering the model, and in section V we conclude with a discussion. 



II. THE POISSON EQUATION IN THE SYSTEM WITH SEC 

Here we reformulate the derivation given in Refi, using the same ideology of separating 
the electrons in core and valence states. The zero level Hamiltonian is the one considered 
in Ref.i, but within a non-orthogonal basis set and, what is more important, here we avoid 
the slave-boson technique in favor of the diagram technique from the atomic limit developed 
in Ref.Si. Let us consider an ion which has n /-electrons in the ground state. Then, only 
the transitions T n —>■ T n ±i will be allowed in the spectrum of excitations while all other 
transitions like T n — > T n ± 2 , T n ± 3 , involving a larger number of electrons, will be strongly 
suppressed by the large energy separation between these states. If the energy of the atomic- 
like transition A 2 = E^ 1 ^ — E^ between an (n + 1) and n-electron state, T and V, of the 
/-ion is much higher than the Fermi energy, Ep, the number of /-electrons in the ion will 
be fixed. Indeed, in this limit this upper "single-electron" level is empty while the lower 
one, even if it forms a band, will be fully filled. In the rare earth elements the populated 
part of the /-spectral density corresponding to the transitions Ai = Ei n ) — E< n -\) is much 
below Ep (except for Ce, and perhaps Sm). It can be also much below the bottom of the 
conduction-electron bands. When an orbital has such a low energy, the mixing interaction, 
as well as overlap between these core-like levels and conduction electrons are negligible^. 
This physical picture exactly corresponds to the type of ab initio calculation where the /- 
electrons are kept in the core. One can use the many-electron functions for the description of 
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the ground state of an ion from DFT-LDA-based calculations too. All electrons in this case 
experience the same potential. As has been shown in Ref.l, this picture is valid when the 
energy of the upper transition, A2 = -E^n+i) ~E n , is much above Fermi energy. However, the 
photo-electron spectroscopy experiments show that even in rare earth elements, for which 
this picture seems to be most appropriate, the level A 2 is sometimes only slightly above 
Therefore, due to mixing interaction and, possibly, hopping, a band with mixed /- 
and conduction electron statesEl can be formed. As discussed in details in Ref.l, this leads 
to shifts of spectral weights from integer values and a violation of the single-electron picture. 
These spectral weights, therefore, must enter the expression for the charge density in the 
Poisson equation. Besides, they control the strength of the mixing and hopping. Let us 
derive the Poisson equation which contains the information about these spectral weights 
and is valid not only for zero temperature. Here, we will not consider the formation of 
phonons and scattering processes which involve them, nor plasmons. 



A. The Hamiltonian in a Non-Orthogonal Basis Set 

An orthogonalization procedure of the wave-functions belonging to different ions leads 
to a coupling of the states. This makes it difficult to separate the strongest single-site in- 
teractions. Therefore, the local strong interactions between /-electrons can most easily be 
taken into account in the non- orthogonal site representation. For this reason we will, to 
some extent, use the technique developed previously! (below referred to as I). The delocal- 
ized electrons are treated within the weak-coupling perturbation theory (WCPT), while the 
localized (or semi-localized) within the strong-coupling theory (SCPT), see I. In order to 
introduce, for the /-operators (and other core electrons), the many-electron representation 
we rewrite the field operator, ipo-(r), in the jL-representation 



W(r) = J dqe- l ^<P jL (r)a jL , (1) 
a iL = [(1 - fo,»)c jL + S L>ll f jlt ). (2) 
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Here, j = Rj is the site, L = (I, mi, s = 1/2, a), I is the orbital moment, mi is its projection 
to the z-axis, s is electron spin and a its projection to the same axis /i indicates localized 
electrons. In Eq.(|2]) we have separated all electrons to two classes: core electrons, fj^, which 
either remain fully localized in solids, or only partly delocalized, and delocalized electrons, 
CjL, which will be described in k-space in regular crystals. Since the essential part of the 
work to be done concerns the localized electrons, it is reasonable to formulate the approach 
in the site representation first. The basis functions 4>jl(?) are in general not orthogonal to 
each other, 



J dr0* L (r)0 i/L ,(r) = O jL # v . 



(3) 



therefore, from {^(r), ^(r')} = S(r — r')5 aa ', we see that 



where Ojl,, L , is the (JL, j'Z/)-matrix element of the inverse of the overlap matrix O. 



(4) 



The full Hamiltonian is 
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C n 
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(5) 
(6) 



where is the infinite constant J drdr'^r — r')v(r — r') J2a' &aa' which arose when we 
transformed iplip2 v (^-i 2)V , 2 , 0i m to ip\ipiv(l, 2)^2^2- We omit this constant below, since it 
does not influence the physics. Let us rewrite the Hamiltonian in the representation using 
the functions 4>jL(?) (which can also be defined in different ways and we shall discuss it 
later). Using the expansion in Eq.([T|) inserted into Eq.(^) gives 



.air <l ;.. 



Z {jL} 



(7) 



Here 



32 l 2 J3L3 



U2L2 



2m 



\hLz 
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Now we assume that the nuclei are in fixed positions and separate the part of the Hamiltonian 
which contains on-site interactions between electrons that are treated as core electrons 



njion \ * njion \ * J \ * r f 



+ -Yv- ■ ■ 



(9) 



The single-site part of the problem and the rest will be treated in different approximations. 
We want to use Hubbard operators that are usually introduced in such a way that they 
diagonalize the single-ion Hamiltonian 



H™\j,r n ) = E jrn \j,r n ), 
xf^ = \j,r n ){j,r r , 



(10) 

(ii) 



Let us now discuss how to construct them. 



B. The Hubbard Operators in Terms of Fermions 

We are not able to diagonalize the full Hamiltonian exactly, unless for certain model cal- 
culations, and we have to use some approximation. In order to ensure that the ground-state 
wave function fulfills Hunds rules, we will follow the technique by IrkhinEl, who translated 
the Racah technique, used in atomic spectroscopy^ for the wave functions, into the oper- 
ator language. Although this issue has been discussed in Refs.ll3i, we find the definition 
used not fully satisfactory and in need of a slight modification. For this reason we discuss 
the definition of the many-electron operators and, correspondingly, the connection between 
the Hubbard operators and the many-electron creation and annihilation operators and the 
modification needed. The creation operator for a group of n equivalent electrons (say, in an 
/-shell) in the many-electron state \T n ) = |0) in Irkhin's definition has the form 

7 



At = —= Y CE" Cf> U /14 • (12) 

V ^,r n _i 

Here are the fractional parentage coefficients, which do not depend on the momentum 

projections (if n < 2, G = 1 and the squared coefficient, (G ! r ™_ 1 ) 2 , measures the fracture of 
the state |r„_x) in |r n )); C r ™_ are the Clebsch-Gordan coefficients, 

° r »-l.f - U L„M|"- 1 ),Zm i U 5„_ 1 M("- 1 ),l/2, ( T' ^ 

where L n , , S n , Mg are the orbital moment, its projection, spin moment and its projec- 
tion for the n-electron configuration |r„). In order to be able to calculate the commutation 
relations between the conduction electrons and the Hubbard operators (H-operators) as well 
as between the H-operators themselves, we have to express them in terms of Fermion oper- 
ators. We have to provide the correct commutation relations for H-operators, belonging to 
the same site 

[X~< r , X Kx ] ± = <5 rA X 7X ± 5 X7 X Ar . (14) 

If we try to define an H-operator in the form of a product of the operators A, i.e. X lV = A^A-p, 
we should get zero if we multiply by X /yTn X AmX for n ^ m. This should be provided by the 
fact, that for the Fermion operators f 2 = (f^) 2 = 0. However, the product Ap n A\ m ^ 0, if 
n < m, because the operators A^ n do not contain information about non-filled orbitals of 
the shell. The operators A Tn do not provide the orthogonality of the states with different 
number of electrons in shell. The recipe suggested in Ref.0 is to define new operators as 
follows 

4„ = 4. IK 1 - 4„ = I^ 1 - K)Ar n - (15) 

The product should run over all non-filled orbitals. Let us consider, for example, the two- 
electron state composed of /-states, |r 2 ) = \L = 5, M K L = 4, S = l,Mg = 0). Then, the 
operator 

4„ = ^ , 3mi,3m2Cl/2,cr 1 ;l/2,cr 2 /3mi;l/2,o-i/3m2;l/2,o-2 



combines three terms (for briefness below we omit indices I = 3 and spin s = 1/2 : 

4 2 = af/sVu + 7/ 2 f T / 2 \ + /f T /U (17) 



with a = J 5/253, 7 = 27/vl5. Now, we have to multiply this by the product of the 
factors (1 —U/j,), where « runs over all empty orbitals. From this a problem is apparent: the 
different terms of the combinations of Clebsch-Gordan coefficients involve different orbitals 
and, therefore, it is impossible to choose a single factor which includes all empty orbitals 
for each term in the sum. Therefore, each term of the sum must be supplemented with its 
own factor. In this particular example the new many-electron operator should be defined as 
follows 
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2 « 



mi7^3 m2^1 m\^2 m 2 ^2 



+ AVi, II ( X ~ "mit) II i 1 - ™rn 2 l) 



a [4 2 (3T,U) + 7^r 2 (2T,2i) + 4 



r 2 (it,3J.)J- 

(18) 



Obviously, in a general case, when we construct the operator Ap 2 , each term in the sum 
should be multiplied by the product of the factors (1 — h ma ) corresponding to this term, 
where the set of {ma} includes only those orbitals which are not included in the product 
of the /-operators in the corresponding term of the Clebsch-Gordan sum. Let us denote 
this product ft, where the bar over IT means that it contains only complementary orbitals. 
In I we used the orbital representation, where each H-operator contains only one term, 
containing n creation /-operators for the n-electron state which is multiplied by the product 
fl. We we will call it elementary operator. In the case of Eq.flT8|) these operators are 

4 2 (3TU)' 4 2 (2t,2j.)j 4 2 (iT3!)' ^ e can ' therefore, make the statement that any operator 
of a state in a central field (i.e. of the Clebsch-Gordan combination type) or in a crystal 
field, can be represented as a sum of elementary operators with the coefficients, which are 
dictated by the symmetry of the surrounding of the ion. Since the construction of the state 
in central field within the Racah technique is recursive, i.e. the n-electron state is composed 
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of (n — l)-electron states and one-electron states, the (n — l)-one is made of the combination 
of (n — 2)-and one-electron states, and so on, this statement needs proof. Let us start with 
a many-electron operator, describing one localized electron in orbital 1. The electron state 
has the form f\Y\ v ^ 1 {l — h u ). In order to construct the state, which contain two electrons 
localized in the states 1 and 2, we have to multiply this operator by f\ from the left-hand 
side. Since /^(l — ^2) — Hi a U extra factors (1 — h u ) are automatically projected out and 
we are back at Eq.(|l2|). Therefore, each step to a higher number electron operator will be 
started again with Eq. (|l2"l) . This means that the factors II should be added in each term of 
the sum in the last step only. Thus, the Hubbard operator can be written in the form 

X^ r = A\A r , (19) 

where each elementary operator entering the sum for the operator contains the projecting 
product flW. The same is valid for Ay and flp . Then, we can represent each /^-operator 
in any place where we meet it, particularly, in the Hamiltonian, in terms of H-operators 

f m EE {f^fXl (20) 

Here, repeated indices are summed over. The H-operators can also be written in the Hubbard 
formX] 7,r] = |R i) n,7)(R i) n+l,r| or in terms of products of Fermion operators as discussed 
above (here |Rj,n, 7) and \Rj,n+ 1,T) are many-electron n- and (n + l)-particle states). 
Using the definition of the many-electron operators and the Hubbard operators in terms 
of Fermion operators one can calculate all commutation relations between the conduction 
electrons operators and the Hubbard operators!!, 

{c jL ,X}} = Ojl j , v fr£4. (21) 

Here, a, b denote the Fermi-like transitions from the n- to (n + l)-electron state (a = 
ct(r n ,r n+ i)), a denotes the inverse transition (a = a(T n+ i,T n )), £ = £(r n ,r^J and e b ^ a 
are the structure constants of the algebra for the Hubbard operators, 
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The summation over repeating indices is implied. A Bose-like transition is denoted as 
and a diagonal Bose-like operator is denoted as hj, hj = Z| ( - r ' r ' ) . 



C. The Field for the Electric Potential and the Poisson Equation 

Let us introduce the Fourier- component of the charge density for the nuclei in the point 
R a with the charge Z a : 

Pn(q) = E^ qRj - (23) 

3 

Then, the interaction between the nuclei is 

*C = \ E I' 2 ''* I = / d ^^f^- (2 4 ) 



The interaction between electrons and nuclei, 

^r = -E/drp e (r)-^ 



Zj e 2 



3 A 



/Arce 2 
dqp e (q)^^Pn(-q), (25) 
Q 



contains the electrons which are localized to the same nucleus, or to other ones, and those 
which are in a mixed state. The interaction of the localized electrons with their own nucleus 
is 

Z,e 2 



-E/ dr^^Kir-R^Mr -nMf^Zl 

j J 

= -/dq^E^f(q)4 (26) 

q 3 



where Sj is the form-factor of the ion, which takes into account the contribution of the 
orbitals v and /i into the transition £ 

0%q) = J dre^<p* v ( r )<P m {r){flf^ (27) 
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and the index j in JAt (r) denotes affiliation of this function to the ion on Rj. This term is 
included to the Hamiltonian for the ion in this point. The term 

-E/dr^-(l-MC(r-Rn) 

4iie 2 
T 2 

describes the interaction of localized electrons belonging to one ion with a nucleus of another 



= - / dq ^T E^e-^Oi^Zie^] (28) 



ion. Then, the operator of an effective density of ions can be written as 

Mq) = E[^-E^(q)^|K qRj (29) 

and the interaction between all nuclei which are screened by their electrons is 

^/dqA(q)^A(-q)-Wo- (30) 

Here, Tin takes into account the terms which are double counted in the first term, since 
there is no interaction of the ion with itself 

Hd = \ E / dq[Z, - £ Of (q)Z}]^[^. 
3 5 q 

-E<?f(q)4l- (3i) 

The terms of this interaction at small q, as well as Hd itself, diverge. This is the standard 
problem of screening. 

All other electrons belong either to the class of the transitions between different ions, or 
to a mixed state between conduction-electron and localized electrons, or to the conduction 
electrons. The operator of the charge density, p r (q), of these remaining electrons can be 
written as follows 

Pr(q) = E OjLfV^jLCfV 

jLj'L' 

+ E ^x /M '(q)(40 a 4^ + ^a'(q)(/5^|^^ 

jLj'L' 

+(! - ^0O iwV (q)(/t) a (/ M 6 AjX» (32) 
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The matrices 0(q) are defined as 

O jL M<l) = J dre^- r 0* L (r)0 J , L ,(r). (33) 

Thus, we obtain the Hamiltonian in the many-electron representation 

H = H + (H coul + T), (34) 

where T is the kinetic energy, and the zero Hamiltonian is 

Ho = Hi + H c = £ E jV h] + J2 h% JIL/ c] L c jlLl . (35) 

Here 7ijf describes the electrons treated as core electrons when the interaction between the 
ions and all other electrons is switched off 



= (j,ri y ^- - y - 

and the function \ j, T) is a many-electron wave-function. The Hamiltonian of the conduction 

electrons, Hq, describes the electrons in the field of the unscreened nuclei. This will be 

supplemented with the self-consistent field from the electrons. The kinetic energy of the 

electrons which belong to one site is taken into account, while the part T of it, which is 

proportional to p r (q), is n °t- We will treat this part of the kinetic energy as a perturbation. 

Thus, the full Coulomb interaction can be written as follows 

1 f Ane 2 

- j dq[p,(q) - p r (q)]— — [p t (q) - p r (q)] 

1 r Aire 2 

= 2J dc ^q^-2-Ar (37) 



The partition function, written in the standard form, is 



Z = Trexp[-/3{H- fiN)] 

= Z ^Tr{exp[-f3Ho]T T exp[- f dr^ int (r)]} 
Z Jo 

= Z ■ {% exp[- f drH mt (r)]) (0) (38) 
Jo 
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Let us write the part H c ° p ul of Hint in the Fourier transformed form, Eg. fl32"|) . Then, we 
introduce the Gaussian functional integral 

1 = J Dif^r) exp[- jf" ^ |dq q 2 ^)^^)} (39) 

under the trace of the partition function. We make the shift 

Attc 

VqM -» VqW + ^^Pq, ( 40 ) 

in this Gaussian integral. This allows us to rewrite the Coulomb interaction in terms of 
interaction of electrons with the random field v 9 q( r )- This shift generates the term (—H C p P ul ) 



which cancels the H c ° p ul in the Hint but adds the term 



ic r@ r 

~2k dT J dq [ ^ (r)/3 - q(r) + ^( r )v 9 -q( r )]- ( 41 ) 

Note that, although the /-orbitals of the same site are orthogonalized, the q 7^ O-components 
of the overlap matrix have non-zero values and, therefore, non-diagonal transitions £ 7^ [T, T] 
enter the expression for q-component of the charge density. Thus, we have to work with the 
following expression for the partition function 

Z f ff 3 f q 2 

7T = J D( Pq{vexp[-J Q dr J dq — ^ q (r)^_ q (r)j 

xexp[-/^ q (r),y;_ q (T)]], (42) 



where 



= -- ln(T T exp[- ^ dr | dq(T q - H 



-H D )])^\ (43) 



T 



^2 

2m 



PM) (44) 

q 

= 2 / dq t( ie< Pq( r ) ~ v)P-q( T ) + Pq( r ) 

x(ie^_ q (r) -m)], (45) 
H° = H C + H%. (46) 
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Here, we have used the fact that the fields y^qfV) commute with any operator. The function 
T is written in the form usually used for the cumulant expansion. 

Since both the mixing interaction and overlap matrices are non-zero, a part of the charge 
is in the mixed c/-states. Besides, the /-subsystem is described in terms of non-linear X- 
operators. For these reasons we cannot describe the full contribution from the term jiN 
in the zero Hamiltonian. The field tpqij) can be interpreted as the field for the electrical 
potential which is generated by Coulomb interaction. If we put e = the system does not 
have charged particles and, therefore, there are no contribution to the partition function from 
this field. At e 7^ in a non-homogeneous system an average charge density is not equal to 
zero, and, therefore, the expectation value (y^qfV))^ is non-zero too. Taking the functional 
derivative of the free energy, we find that the saddle-point approximation generates the 
Poisson equation for this field 

5F _ 5(-T\nZ) 

= -^Wr)> v + |(p q (r)>^ = 0. (47) 

The static part of the field of the electrical potential $(r) is connected with the field (</?) as 
follows 

$(r)= lim i( V (r,r))^i( V (r,0)) r (48) 

T— » — 

Thus, we have the Poisson equation for the electric field $(r), 

V 2 $(r) = -47r(p q (0)>„, (49) 

where the charge density is the difference between the densities of the delocalized electrons 
and "soft" ion (where the polarization of the ion and excitations are allowed). 

III. THE SELF-CONSISTENT HUBBARD- ANDERSON MODEL 

Since the charge density, (/^(O))^, contains expectation values (Cj L Cy L ,), (Cj L Xj,), 
{Xj,c jL ), (XjXj,), we need expressions for them. Since the field </?(r, r) contains the av- 
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erage static field and deviations of it 

(P(t,t) = j$(t) + 6<p{t,t), (50) 

it is reasonable to start with the approximation iip^r)^ — > «((^>q(0)) ¥ , — > $(q), and to neglect 
the fluctuations of this field. Then, we obtain the following saddle-point Hamiltonian 

H = jdq g 2 $(q)$(-q) 

+ J2( h JL,j'L' ~ fi OjL,j>L>)c] L C jlL , 

+ E(^r ) -n/i)^ r , (51) 
Hint = T + C -J dq [$(q)p_ q + p q $(-q)] - (52) 

Here 

h jL , j/L , = (jL\^- + $(r)\fL') (53) 

is the frequency matrix of the conduction electrons in the self-consistent field <£>(r). The 
Hamiltonian h for them can be obtained if we make a transformation to the orthogonal 
variables a, using the Cholesky decomposition for the overlap matrix O: 

c f (h - 0)c = ctfZtZ^hZ-^Z-^ ZZ]c = 
(c f Z) [h-pl] (Zc) = a f [h-/il]a. (54) 

Let us now insert into Hint the expression for p q in the jL-representation. Taking into 
account that 

| / dq S(q)0 iW ,(q) = | / dr 0* L (r)$(r)0 ja ,(r) 

= ^LJ'L' (55) 

are matrix elements of the self-consistent Coulomb field, we find that Ti actually gives the 
periodic Hubbard-Anderson Hamiltonian 
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H = y dq g 2 $(q)$(-q) 

+J2(h LDA -»o) jL!fL ,c] LCj , L , 
+ E[W(fe)^ [r ' rl ^ + ^.c.] 

+E(i-%0(|^ + ^-^W,' 

x(f})*(f„,) h X}X} 

-H D . (56) 

Thus, within this approximation, the single-ion energies E^ are shifted by the self-consistent 
field of interaction with other ions, delocalized electrons and the localized ones, but belonging 
to other sites: E^* = + (r|$(r)|r). However, the self-consistent field should be found 
from the Poisson equation and, therefore, it contains only the Hartree part of the interaction 
between the collective quasi-particles. Nevertheless, this problem does not coincide with the 
single-particle Hartree approximation, since it contains additional information about the 
structure of the many-electron states |T) of the ion. Due to the non-linearity of the problem 
the approximations for the Green functions (GFs) can be used in a different form from 
the standard single-particle problem in the Hartree approximation. This does not only 
lead to different expressions for the charge density, but makes it possible to improve the 
approximations for the matrix elements. Thus, in spite of the fact that the static field $(r) 
fulfills the Poisson equation, in this approximation VjLyv does not necessarily correspond to 
the potential for the electrons in the Hartree approximation. Particularly, as seen from the 
solution of the self-consistent Anderson model, presented in Ref.0, the potential depends 
implicitly on the many-body population numbers of the ion states already in the lowest 
approximation. 

The exchange contribution appears in first order between electrons via the fluctuation of 
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the field <p, oc (T5ip(r)5ip(r')). This study we leave for the future, however, it is interesting 
to note that this exchange involves also the contributions from fluctuations caused by the 
intra-ion transitions. 

The charge density for the Poisson equation can be found from the GFs. It is clear that if 
it is possible to approximate this potential by a spherically symmetrical one, V} w y = S^'Vj, 
then (f^fnY = %[r,r] for the orbitals, /i, occupied in the state T. The averages (c^c), (c^X), 
(X'c), (X'X) needed for finding (p(r)), can be found from the corresponding GFs. 



IV. MIXING INTERACTION 

Let us now rewrite the problem in a form close to the standard periodical Anderson 
model. This allows us to discuss different possible definitions for the matrix element of 
the mixing interaction. Usually the Hamiltonian of the Anderson model is written in the 
following form 

H (An d) = J2elclc ka + E[^;e ik - R ^L/ JM + H.c.) 

kcr kcr 

+ Um^nsmf jmf j/Jbf j^f 3H4, ^2 e t J 'fj^fj^' (^) 

Here is the spectrum of the conduction electrons, V£ is the matrix element of mixing 
interaction, 

2 

= / dr ^W[|^ + V(r))cp„(r - R,), (58) 

and U^n^m is the matrix elements of the Coulomb intra-shell interactions. For the mixing 
no problem arises in the case of the impurity Anderson model, since the potential in this 
case is the difference between the periodical potential for the conduction electrons and the 
potential of the impurity, although the local on-site term is equal to zero due to symmetry 
reasons, overlap with the orbitals of neighboring ions gives non-zero contribution (see the 
paper by Anderson0). In the case of the periodic Anderson model (PAM) this form of the 
matrix element suggests that either the potential has different symmetry from the Bloch 



wave functions, or the functions ^{r) and </v( r — ^j) are n °t orthogonal. It is not clear 
how to fulfill the first assumption in the case of elemental metals (like Ce metal) since we 
are dealing with a periodic system. In the second case 



Therefore, in order to have a non-zero mixing, one has to work in terms of an non-orthogonal 
basis set, but in this case a) there is a contribution from the chemical potential to the 
partition function which affects the mixing, and b) the non-orthogonality causes non-zero 
anticommutation relations between the /- and c-operators. This is never taken into account 
in model calculations. 

Let us consider the consequences from our formulation. Since the mixing matrix element 
has a single-electron form, let us make in our saddle-point Hamiltonian a transformation that 
diagonalizes the conduction electrons. First we have to rewrite the Hamiltonian, Eq.(|56|) in 
k-space (for the formulas to be transparent, we will write them for the case of one atom in 
the elementary cell, then, the single-ion matrix elements and energies do not depend on the 
ion index). Using the definition: 



/dr0r(r)[^ + nr)]^(r-R,) 




(59) 




(60) 



k 



A' 



we find 




+ J2[hL,L>(k)-nO L>u (k)]c } 



; kL C kL' 



+ J2(E^ - n^)h] + U^Aflf,') mT ' T] Zl + h.c] 

+E(f^ + v- 

+ E(f^ +V- ^W^(4W 

+E(i - W(£ + v- 



—TCd- 



(61) 
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Now we have to rewrite this in terms of the operators a k7 = Z^ L c kL , a kl = c kL Z k 
easy to see that they are orthogonal to each other 

«k 7 «ky + «k 7 '«k 7 = z l L ( c ^L c iv + C kL' C kL)^k V 

= zi L (o- 1 ) LL ,z^' = zl L {z^{z-y^z^' 

c f (h - n 0)c = c^Z^hZ-^Z-// ZZ]c = 
(c t Z)[h- / uI](Zc) = a^[h-fil]a. 

Thus, we get 

H = Jdq g 2 $(q)$(-q) 

+ E[(4" 1 ) 7L ^(k)(^ k - 1 ) i7 ' - /i5 7 y]4 7 « k y 

+ - n/i)^ r + EtW^O^'^l + h.c] 

+ UZ- 1 r L (^ + V- ^0) jW e^(f^lXJ, 

+E(f^ + ^ - ^^(zir 1 )^^)^-*^^ 

+E(i - M(f^ + ^ - vo^ifinf^x}*!' 

—T~Cd- 

At last, diagonalizing the conduction electron Hamiltonian Z~ 1 hZ~ 1 , we have 

H = J dq g 2 $(q)$(-q) + 5>kA - A*]4^A 
+ E(4 n) - n/i)^ r + E[^,'(4/,0^ [r ' r] ^| + h.c] 

+E^K k M A ^ 

+ Ee- Jk - R ^X|a k7 

+ E(i - ^)(|^ + ^ - MwAfinfsn!*} 

—7~Ld, 

where the band energy and mixing matrix elements are 
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(66) 



VT(k) = #;\k)(Z-y L (i- +V- /^) iW (4' 



2m 



i(k)(/ M O a , 

2 



^"( k ) = (£r + v ~ »o)i»,?L{z?)^e^){fl) 



2m 
(4)^?(k), 



(67) 



(68) 



where {■$} diagonalizes Z _1 hZ _1 , and v£ is the one-electron hybridization parameter Thus, 
the matrix element of the mixing interaction has to be found self-consistently together with 
the charge density (which in turn depends on the particular approximation in which the 
PAM is solved) and can be represented in the form of a sum over all localized orbitals of 
products of the matrix element of the single-particle potential (on the conduction electron 
Bloch function 4>k\{ r ) and the localized orbital x^i r )) an d the factor {f^') a , which reflects 
the contribution of this orbital into Fermion-like transition. IrkhinEl has performed the 
calculation of this factor for the 4/-elements making use of the Racah technique leaving, 
the single-particle matrix element undefined. Putting a = a(T n , r n _i) and \i = (I, m, a) and 
using the result of the calculation of Irkhin, we can write these coefficients as follows 



(flT = <r B |/L,|r„-i> 



n[T n ][T 



Sn 



Jn 



n+1 



} G T n ! 

J- 71 — 1 



(69) 



Sn— 1 L n —\ J n -\ 

1/2 I j 

where 7 = (Zm/cr) are the one-electron quantum numbers, [a] = 2a + 1 and G^ n is the 
parentage Racah coefficients and j = I + 1/2. The coefficients (/^) a are different for the two 
channels, j = 5/2 and j = 7/2. The square of them are given in Table. |I[ 
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TABLES 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


L/5/2] 2 


1 


2.25 


2.3 


1.8 


0.9 


0.15 





1.6 


2.3 


4.1 


1.7 


0.9 


0.3 





[f7/2? 





0.25 


0.6 


0.8 


0.63 





0.9 


0.5 


1.3 


4.1 


3.0 


3.3 


2.8 


1 



TABLE I. Coefficients for the mixing interaction for many-body states. For the /-series, taken 



fromlll. 
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However, there are more considerations to the problem, since the Bloch orbitals are not 
orthogonal to the localized orbitals. In the equation of motion for the operators appear 
combinations involving the overlap matrices. For the Hamiltonian which includes hopping 
and mixing interactions (the periodic Hubbard- Anderson model) the diagram technique and 
different approximations for the Green functions (GF) are given in I, while the full self- 
consistent solution for rare earths in the simplest possible approximation is presented in 
Reffi 

V. DISCUSSION 

Many different suggestions exists in the literature on how the parameters of the Ander- 
son Hamiltonian (periodical or impurity) should be calculated from an ab initio approach 
and a consensus is not yet achieved. This motivated us to make an attempt to derive the 
parameters from a total Hamiltonian. Since the Anderson model is usually used for the 
description of the cases when strong electron correlations are well developed, it is reason- 
able to consider strong intra-atomic Coulomb interactions first. For this reason we have 
performed a derivation of the self-consistent PAM in four steps: i) first we separate, in 
the Hamiltonian, the strong intra-atomic interactions and approximately diagonalized them 
with the help of many-electron functions describing different ion terms; ii) we expressed all 
operators, describing /- and core electrons in terms of Hubbard operators; iii) we performed 
the Hubbard-Stratanovich decoupling of the Coulomb interaction and finally iv) we found 
the equation for the electric potential in a saddle-point approximation. This lead us to the 
effective Hamiltonian, which coincide with the generalized PAM. 

All three operations before making the saddle-point approximation are exact. However, 
the saddle-point approximation neglects the contributions from the exchange interaction for 
the delocalized particles (excitations) and Coulomb screening effects, which appear only in 
the next orders with respect to fluctuations of this field near its saddle-point value. Thus, we 
may conclude that from this point of view this model is quite rough. On the other hand, the 
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way suggested here has an obvious advantage compared to the PAM with non-self-consistent 
parameters since it, at least, takes care about perturbations of the local charge density which 
may introduce quite large changes in energy. 
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